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ABSTRACT 


We report on a short XMM-Newton observation of the radio-quiet Narrow 
Line Seyfert 1 PG 1402+261. The EPIC X-ray spectrum of PG 1402+261 shows 
a strong excess of counts between 6—9 keV in the rest frame. This feature can 
be modeled by an unusually strong (equivalent width 2 keV) and very broad 
( F WHM velocity of 110000 km s 1 2 3 4 ) iron 14-shell emission line. The line centroid 
energy at 7.3 keV appears blue-shifted with respect to the iron Ka emission band 
between 6.4 — 6.97 keV, whilst the blue-wing of the line extends to 9 keV in the 
quasar rest frame. The line profile can be fitted by reflection from the inner 
accretion disk, but an inclination angle of > 60° is required to model the extreme 
blue-wing of the line. Furthermore the extreme strength of the line requires a 
geometry whereby the hard X-ray emission from PG 1402+261 above 2 keV is 
dominated by the pure-reflection component from the disk, whilst little or none 
of the direct hard power-law is observed. Alternatively the spectrum above 2 
keV may instead be explained by an ionized absorber, if the column density is 
sufficiently high (Ah > 3 x 10 23 cm -2 ) and if the matter is ionized enough to 
produce a deep (r ~ 1) iron K-shell absorption edge at 9 keV. This absorber 
could originate in a large column density, high velocity outflow, perhaps similar 
to those which appear to be observed in several other high accretion rate AGN. 
Further observations, especially at higher spectral resolution, are required to 
distinguish between the accretion disk reflection or outflow scenarios. 


Subject headings: galaxies: active — Seyferts: individual: PG 1402+261 
rays: galaxies 
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1. Introduction 

PG 1402+261 is a type-I radio-quiet quasar at z=0.164 (Schmidt & Green 1983). It emits 
narrow optical permitted lines with FHWM (H/3)=1910kms -1 and can also be classified as 
a Narrow Line Seyfert 1 (NLSl). The mass of the supermassive black hole in this quasar 
has been estimated to be 2xl0 7 M© and the bolometric luminosity is about 1.3xl0 45 ergs -1 
(Woo & Urry 2002). This leads to a high accretion rate relative to Eddington of about 50% 
(Porquet et al. 2004). PG 1402+261 was previously observed in X-rays with Beppo-SAX by 
Mineo et al. (2000). They found that the data were satisfactorily represented by a broken 
power law continuum of T = 2.59 ± 0.10 and T = 1.52 ± 0.30 for the soft and hard X-ray 
bands, respectively, although there was some evidence for an excess of counts above 7 keV in 
the MECS spectrum. A preliminary analysis of a recent (2002) XMM-Newton observation 
of PG 1402+261, reported by Porquet et al. (2004) amongst a sample of PG quasars, shows 
that the spectrum cannot be represented by a simple broken power-law continuum alone, 
but exhibits a large deviation between 6-9 keV in the quasar rest frame, in the form 
of an apparent broad and blue-shifted iron K emission line. This observation adds to the 
number of AGN (and NLSls in particular) that appear to exhibit unusually strong emission 
or absorption features above 7 keV in the iron K-shell band; e.g. 1H 0707-495 (Boiler et al. 
2002), IRAS 13224-3809 (Boiler et al. 2003), PDS 456 (Reeves, O’Brien & Ward 2003) and 
PG 1211+143 (Pounds et al. 2003). These iron K band features have either been interpreted 
as the signature of enhanced reflection from the inner accretion disk, e.g. 1H 0707-495, 
(Fabian et al. 2002), but in some cases may also be interpreted as evidence for an additional 
high column density absorber, either from matter partially covering the X-ray source (Boiler 
et al. 2002) or even from a high velocity (~ 0.1c) outflow (Pounds et al. 2003; Reeves, O’Brien 
& Ward 2003). 

Here we present a spectral analysis of the XMM-Newton observation of PG 1402+261, 
modeling the strong iron K-band emission (or absorption) feature. We investigate the possi- 
ble physical origins for this feature; either in the form of an extreme broadened iron emission 
line produced via X-ray reflection from accretion disk or alternatively a very high column 
density ionized absorber. Note that all fit parameters are given in the quasar rest frame, 
with values of Ho=75kms -1 Mpc -1 , and q 0 =0.5 assumed throughout. Errors are quoted at 
90% confidence (e.g. A% 2 = 2.7, for 1 parameter of interest). 


2. The XMM-Newton observation 

PG 1402+261 was observed by XMM-Newton on 2002 February 27 (OBSID: 0109081001), 
with an exposure time of only 9.1 ks. Data were taken with the EPIC pn and MOS CCDs 
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(Striider et al. 2001; Turner et al. 2001), in Large Window Mode with the thin filter. The 
data was reduced using version 5.4.1 of the XMM-SAS software. Data were selected using 
event patterns 0-4 and 0-12 for pn and MOS respectively, whilst only good X-ray events 
(with ‘FLAG=0’) were included. The source spectra were extracted from a circular source 
region of ~ 40" radius, whilst background spectra were extracted from a box of 1.85' radius 
around the source, but excluding the circular source region. Response matrices and ancil- 
lary response files were generated using the SAS tasks rmfgen and arfgen respectively. 
EPIC-pn lightcurves were also extracted over two energy bands, however no dramatic vari- 
ability was observed, therefore the spectral analysis is performed on the whole XMM-Newton 
exposure.. 


3. X-ray spectral analysis 

In all subsequent fits, the column density is fixed to the Galactic value, i.e, 1.47 x 
10 20 cm -2 . Both the EPIC-pn and the co-added MOS spectra are shown in Figure 1, with 
the data/model ratio residuals plotted against a power-law continuum of photon index 
T = 2.2 fitted in the 2 — 5 keV band. A strong soft excess is present below 2 keV and 
the continuum shape can be parameterized by a simple broken power-law, with a soft in- 
dex of T = 2.84 ± 0.04 and a hard index of T = 2.17 ± 0.12 above a break energy of 
E = 1.8 ± 0.2 keV. A hard excess is also clearly observed above 5 keV in both the pn and 
MOS detectors. For this broken-powerlaw continuum, the unabsorbed 0.3-10 keV (2-10 keV) 
band flux is 7.3xl0 -12 ergcm~ 2 s -1 (1.8xl0 -12 erg cm -2 s -1 ), with a corresponding luminos- 
ity of 4.4 xlO 44 ergs -1 (l.OxlO^ergs -1 ). About 77% of the flux is emitted below 2keV, 
consistent with the strong soft X-ray excess. Note that at this flux level, photon pile-up of 
the source X-ray spectrum is negligible. 


As the pn spectrum contains significantly more counts than the MOS above 6 keV in the 
iron K band, we now concentrate on the pn spectral analysis. Figure 2 shows the data/model 
ratio to the broken-power-law continuum, plotted from 3-11 keV in the quasar rest frame. 
Notice that the majority of the flux of the iron K band feature is observed above 7 keV, i.e. 
above the rest energy for H-like iron, whilst a sharp drop in the spectrum is observed at 9 
keV in the rest frame. A single broad Gaussian emission line models the hard excess rather 
well, with a centroid energy of E = 7.3lo;s keV, a line width of a = keV (FWHM 

velocity of 1.1 x 10 5 km s -1 ) and a very high equivalent width of ~ 2 keV (Table 1, fit 
1). The fit statistic is acceptable, with x 2 /dof = 390.5/377, with the feature detected at 
> 99.99% significance according to an F-test. 
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3.1. Comparison with a previous ASCA observation. 

PG 1402+261 was also previously observed by ASCA on January 12, 2000. The ASCA 
data are composed of the SISO (exposure time ~ 43.2 ks), SISl (~ 37.7ks), GIS2 and GIS3 
(each ~ 57.3 ks) data. A fit above 1 keV using all four detectors, with a single power-law 
(r = 2.3 ± 0.2) plus a broad emission Gaussian line, gives a good representation of the data 
(x 2 /d.o.f=404. 2/361) and the line is required at a confidence level of >99.95%. The energy 
of the line is 7.2lo;gkeV and its width is high with cr=0.8+o;4 keV (Table 1, fit 2). The EW 
of ~ 1 keV, is lower than the one found during the XMM-Newton observation, but is com- 
patible within the errors. During the ASCA observation the 2-10 keV unabsorbed flux was 
2.5xl0 -12 ergcm _2 s _1 (SISO), i.e. about 40% higher than the XMM-Newton observation. 


3.2. A broad disk emission line? 

Given the large width of the emission line in the XMM-Newton (and ASCA) fits, we 
attempted to fit the iron line profile with emission from the inner accretion disk around a 
Kerr black hole, using the xspec laor model (Laor 1991). We assumed an inner radius of 
1.2 R g ( R g = GM/c 2 , i.e. one gravitational radius) and an outer radius of 400 R g . Initially we 
fixed the disk inclination angle at 30°, reasonable for a type I AGN. We obtained a very high 
line energy at 8.5l°;2 keV, required to fit the blue-wing of the line between 8.5 — 9 keV, whilst 
the line equivalent width is extremely large at ~ 4 keV (Table 1, fit 3). This line energy 
is much larger than the maximum energy corresponding to H-like iron (6.97 keV), therefore 
this could indicate a significant outflow of the medium producing the line, i.e. 0.3-0. 4 c 
(90 000-120 000 km s _1 ), for H-like or neutral iron respectively. As an alternative model, we 
free the disk inclination angle, but instead fix the rest energy of the iron line emission at 
6.97 keV, i.e. corresponding to H-like iron. This then requires a large disk inclination of 
67^ degrees in order to fit the blue-wing of the line, although the fit is equally acceptable 
(Table 1, fit 4). Thus either a substantial outflow of the line emitting matter is required, or 
the disk has to be highly inclined to the line of sight. 

Given the evidence for a very strong broad line in PG 1402+261, we tested whether the 
spectrum could be fitted with an ionized disk reflection model. We use the model xion from 
Nayakshin, Kazanas, & Kallman (2000) in the simplest lamppost configuration to model 
the emission above 2 keV, together with a soft featureless power-law component to model 
the steep soft X-ray excess in PG 1402+261. We fixed the height of the source above the 
disk to 10 R g , the inner and outer disk radii to AR g (the minimum the model allows) and 
100-Rg respectively. The accretion rate was fixed to 50% of Eddington, whilst the ratio of 
X-ray to disc (UV) luminosity was set to 10% (Porquet et al. 2004). We find the reflection 
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model fits the spectrum well with r so ft=3.1llo;og and r hard =1.70lo;44 for the continuum 
photon indices and a fit statistic of x 2 /do/ = 397/379. The unfolded spectral fit is shown in 
Figure 3. Although the model produces mainly He and H-like Fe emission (at 6.7 keV and 
6.97 keV respectively), a high disk inclination of greater than >60° is required to model the 
extreme blue-shift of the line. In addition, a high iron abundance of 5 x solar (the maximum 
permitted value) is required, whilst in order to model the extreme strength of the line, the 
hard emission component is required to be reflection dominated, i.e. the illuminating hard 
power-law component cannot be directly observed. 


3.3. Highly Ionized Absorption in PG 1402+261? 

As an alternative interpretation, we have investigated whether the strong iron K feature 
can be modeled with a highly ionized absorber. One can fit the spectrum reasonably well 
with a continuum modeled by two power-laws (r so ft = 3.0 ± 0.2 and Thard = 0.9 ± 0.5), 
although the hard power-law has to be very flat to model the hard excess above 5 keV. The 
addition of an absorption edge at 9.0 ± 0.2 keV, with a depth of r = 1.1 ± 0.5, models the 
sharp drop at 9 keV in the rest frame (Figure 2), with a fit statistic of x 2 /dof = 395/378. 
The energy of the K-shell edge then either corresponds to H-like iron (at 9.27 keV) or He-like 
iron (at 8.75 keV). 

We then proceed to fit the data with the xspec warm absorber model absori (Done et 
al. 1992). The model is applied so that the soft X-ray power-law (which may for instance 
originate directly from the disk) is unattenuated, but so the hard X-ray power-law is modified 
by the absorber (i.e. a partial coverer). In order to model the depth of the edge a large 
column density of Nh ~ 2 x 10 24 cm -2 is required. Formally, the 90from this model is 
log £ > 3.5 erg cm s -1 , with an outflow velocity of < 20000 km s _1 , whilst the column 
density must exceed Ah > 3 x 10 23 cm -2 (assuming 5x solar abundance). The best fit 
model is shown in Figure 4, whilst the fit obtained is equally as good as the disk reflection 
model ( x 2 /dof = 396/378). Inclusion of the warm absorber results in there now only being 
a requirement for a narrow iron Ka emission line, with an equivalent width of ~ 600 eV. The 
line centroid energy is 7.4 ±0.1 keV, which indicates that there may still be some blue-shift 
of this component; e.g., ~ 18000 km s _1 compared to the rest-frame energy of H-like Fe at 
6.97 keV. 
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4. Discussion and Conclusions 

The X-ray spectrum of PG 1402+261 shows a strong excess in the iron K-shell band 
between 6 — 9 keV. If this is due to a relativistic disk line, then its parameters are rather 
extreme. The line equivalent width is ~ 2 keV, even stronger than the broad iron line 
observed in MCG -6-30-15 (Tanaka et al. 1995; Wilms et al. 2001; Vaughan & Fabian 2004). 
Such a strong iron Ka emission line can only be produced if the hard X-ray emission in 
PG 1402+261 is dominated by the disk reflection component, with little of the direct hard 
power-law being observed. Furthermore the rest-energy of the line, peaking above 7 keV 
and extending up to 9 keV, implies that the accretion disk is likely to be observed at large 
inclination angles, i.e. > 60°. Such a high inclination would be appear unusual for a type I 
AGN in the context of AGN unification schemes (Antonucci 1993). 

It is possible that large density perturbations in the disk could explain the strong re- 
flection dominated emission above 2 keV. Variations of up to xlOO in density have predicted 
from simulations of highly turbulent disks (Turner, Stone & Sano 2002). Thus a natural 
geometry could occur whereby the direct hard X-ray emission (e.g. from compact coronal 
flares) is obscured by dense, turbulent clouds within the disk surface layers, the same clouds 
then reflect the X-rays into the line of sight (Guilbert & Rees 1988). A similar geometry 
was applied by Fabian et al. (2002) to the Narrow Line Seyfert 1, 1H 0707-495, which shows 
an unusually strong iron K edge at 7.1 keV. Fabian et al. (2002) propose that the hard 
X-rays could be emitted within low density cavities between dense rings of matter on the 
disk surface; the ring surfaces then reflect the X-rays into the observers line of sight, whilst 
the direct emission is obscured. In PG 1402+261, the direct X-rays are more likely to be 
obscured if the disk was viewed at large (side-on) inclination angles, as is inferred from our 
accretion disk fits. The soft X-ray (and EUV) excess emission may not be obscured however, 
if this is the Comptonized thermal emission from the extended inner disk region (Malkan & 
Sargent 1982; Czerny & Elvis 1987). 

One other possibility is that the disk reflection is enhanced by the effects of gravita- 
tional light bending of the hard X-ray emission towards the disc, i.e. returning radiation 
(Cunningham 1975). This effect will be strongest when the X-rays are concentrated within 
2 — 3 R g around the spin axis of a maximally rotating black hole (in the so-called ’’lamp- 
post” configuration). Indeed it has been suggested that the strong, redshifted iron line in 
MCG -6-30-15 may be boosted by this process (Martocchia, Matt & Karas 2002; Fabian 
& Vaughan 2003). However in the case of PG 1402+261, it may be difficult to explain the 
observed blueshift of the iron line in this scenario, as one would expect the emergent line 
profile to be gravitationally redshifted below 6 keV. 

Alternatively, the spectrum of PG 1402+261 can also be modeled by a high column 
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density, high ionization absorber, where the sharp drop at 9 keV in the spectrum can be 
fitted by a high ionization (He or H-like) K-shell edge of iron. Both the column density and 
ionization of the absorber have to be rather extreme (Nh > 3 x 10 23 cm -2 and log £ > 
3.5 erg cm -1 s -1 respectively), to produce a r ~ 1 iron K-shell edge at 9 keV. However 
other high column density iron K-shell absorbers, with similar absorber parameters, have 
recently been observed in XMM-Newton and Chandra observations of several high luminosity 
AGN; in particular PDS 456 (Reeves, O’Brien & Ward 2003), PG 1211+143 (Pounds et al. 
2003), PG 1115+080 (Chartas, Brandt & Gallagher 2003) and APM 08279+5255 (Chartas 
et al. 2002; Hasinger, Schartel, & Komossa 2002). The outflow velocities inferred in these 
absorbers were generally found to be extreme, of up to 0.1c. One possibility is that the 
absorber arises from the innermost part of an high velocity accretion disk outflow, within a 
few (10 — 100R S ) radii of the black hole. Indeed as noted by King h Pounds (2003), high 
velocity, large column density disk winds may be a natural consequence of accretion at a high 
fraction of the Eddington rate; in PG 1402+261 the accretion rate is ~ 50% of Eddington 
(Porquet et al. 2004). 

Determining which model (reflection vs absorber) is the most plausible explanation for 
the iron K profile in PG 1402+261 is not possible with the current short XMM-Newton 
observation. However a high column density absorber may impart discrete spectral features 
that could be detected with a longer observation (with XMM-Newton) or with much improved 
spectral resolution (e.g. with the calorimeter-based XRS detector on Astro-E2, due for launch 
in 2005). For instance, resonant K-shell absorption lines from Fe xxv or Fe xxvi may be 
produced, whilst strong absorption lines may also be detected from other elements (e.g. Mg, 
Si and S), making it possible to constrain the kinematics of the outflow. Indeed strong, 
blue-shifted (by up to 0.1c) absorption line were detected in the XMM-Newton observations 
of PG 1211+143 (Pounds et al. 2003), whilst resonant absorption lines may even affect the 
iron line profiles of some nearby Seyfert 1 galaxies; e.g., NGC 3783 (Reeves et al. 2004). 
Alternatively if the spectrum of PG 1402+261 can be explained by disk reflection, then it 
would be highly desirable to determine how the strong iron line changes with flux, in a much 
longer observation. 
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Fig. 1. — XMM-Newton EPIC-pn (black) and EPIC-MOS (grey) spectrum of PG 1402+261. 
Energy is plotted in the observers frame. The top panel shows the datapoints (crosses) and 
the model fitted to the data (solid line), folded through the detector response. The bottom 
panel shows the ratio of the data to a power-law continuum model, fitted between 2-5 keV, 
of index T — 2.2. A soft X-ray excess is clearly seen below 1.5 keV, whilst an excess of counts 
is observed in both detectors above 5 keV in the iron K-shell band. 




Data/Model Ratio 
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Rest Energy (keV) 

Fig. 2. — The iron line profile of PG 1402+261, showing the ratio of the EPIC-pn data to 
a broken power-law continuum fit, as described in the text. Energy is plotted in the quasar 
rest frame. An excess of counts, which can be modeled as a broad, but blue-shifted, iron 
K-shell line is observed between 6-9 keV, whilst a sharp drop is observed in the data at 9 
keV. 
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Observed Energy (keV) 

Fig. 3. — An ionized disc reflection model (XION), fitted to the spectrum of PG 1402+261. 
The crosses show the (unfolded) data-points, the solid line the total model emission, the 
dotted line the soft X-ray power-law and the dashed line the ionized reflection component. 
In order to model the extreme blue-wing of the line, the disc must be highly inclined at 
> 60°. The reflection component must also dominate the emission in the iron K-band (i.e. 
the illuminating hard X-ray power-law is not directly observed) , in order for the very strong 
iron line to be observed. 
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Observed Energy (keV) 

Fig. 4. — An ionized partial covering absorber model fit to the spectrum of PG 1402+261. 
The dotted line shows the unabsorbed soft X-ray power-law, whilst the dashed line shows 
the hard X-ray power-law, modified by the ionized absorber. A high column density of 
> 3 x 10 23 cm -2 is required to model the deep absorption edge present in the spectrum of 
PG 1402+261 at 9 keV. 
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Table 1. Table of iron line spectral fits with relativistic line profiles. The underlying 

continuum is a broken power-law. 


Fit 


CT b or <fox 0 A 

EW* 

X 2 /dof 

Gaussian 

7 0+0.4 
* '^—0,5 

i O+0.6 

i - Z - 0.4 

o n+ L3 
z - u -0.8 

390.5/377 

Gaussian^ 

7 O+0.4 
' * — 0.6 

0 
bo 

1 + 
oo 
if*- + 

i * U -0.5 

404.2/361 

laor 

q r+0.3 

O-3-0.2 

o O+0.6 

4 4+2.0 
+•+-1.7 

389.6/377 



30° (f) 



LAOR 

6.97 (f) 

9 O+0.6 

z ' a -0.3 

67 +5o 

9 i +0.8 

A ' l -0.7 

388.2/377 


a Rest frame energy of line in units of keV. 
b Velocity width (a) in units of keV. 
c Emissivity power-law of the accretion disk, 
inclination of the accretion disk in degrees. 
e Fit performed with ASCA data 
f Parameter value is fixed. 


